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We derive general properties of anomalous diffusion and nonexponential relaxation from the theory 
of tempered a-stable processes. Its most important application is to overcome the infinite-moment 
difficulty for the a-stable random operational time r. The tempering results in the existence of 
all moments of r. The subordination by the inverse tempered a-stable process provides diffusion 
(relaxation) that occupies an intermediate place between subdiffusion (Cole-Cole law) and normal 
diffusion (exponential law). Here we obtain explicitly the Fokker-Planck equation, the mean square 
displacement and the relaxation function. This model includes subdiffusion as a particular case. 
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Introduction. — Many studies have been reported on 
the phenomenon of subdiffusion which is typically ob- 
served when, due to dominating influence of traps (see 
[lrHI and reference therein) , the waiting times of random 
walks become a-stable, a < 1, with an infinite mean. 
However, this picture is only an idealization of the physi- 
cal world. In reality the time of trap life can be restricted. 
It should be taken into account that the traps can be lo- 
cated in some spatial regions from which a walker may 
easily escape. Consequently, in a more general represen- 
tation, the random walks start as subdiffusion, but their 
characteristics become very similar to those of normal 
diffusion at large time scales. One of such clear examples 
is a random motion of bright points (BP) associated with 
magnetic fields at the solar photosphere. The BP trans- 
port in the intergranular lanes with times less than 20 
minutes has a subdiffusive character, but the analysis of 
the BP motion supports the normal diffusion behavior for 
larger times. The experimental result is reported in [4|. 
The present paper is just devoted to such a problem. For 
this purpose we are going to apply the tempered a-stable 
processes for the description of diffusion and relaxation. 
In comparison with the purely a-stable process such a 
process has finite moments, but it saves some important 
rudiments of the stable process too [5-7]. Therefore, if 
its inverse process is taken as a subordinator, it provides 
then a diffusive picture occupying an intermediate place 
between subdiffusion and normal diffusion. 

Tempered a-stable process and its inverse. — The 
model of subdiffusion is based on a flexible Montroll- 
Weiss idea on continuous time random walks (CTRW) 
8]. Briefly, the representation of anomalous diffusion 
by means of the CTRW methodology is the following 



(for more details see, for example, 0, Gil)- Consider a 
sequence Ti, i = 1,2,... of non-negative, independent, 
identically distributed (iid) random variables which rep- 
resent waiting-time intervals between subsequent jumps 
of a walker. The random time interval of n jumps in 
space equals T(n) = £™ =1 T, with T(0) = 0. The ran- 
dom number N t of jumps, performed by the walker up to 
time t > 0, is determined by the largest index n for which 
the sum of n interjump time intervals does not exceed 
the observation time t, namely Nt = max{rt : T{n) < t}. 
The position of the walker after Nt jumps becomes then 



R{N t ) = J2Ri 



(1) 



i=i 



where Ri are iid variables giving both the length and the 
direction of the z-th jump. The process ([!} is just known 
as the CTRW. 

If the time intervals Ti belong to the domain of at- 
traction of a completely asymmetric a-stable distribu- 
tion with the index < a < 1, the generalization of 
the central limit theorem yields the continuous limit 
a _1 / Q r([ar]) U T as a — » oo, where U T is a strictly 
increasing a-stable Levy process, a > parameter, [x] 

denotes the integer part of x and "A" means "tends 
in distribution" . Similarly, let the jumps Ri belong 
to the domain of attraction of a 7-stable distribution 

5 7 ,^(x),0 < 7 < 2, \p\ < 1 so that a - lh R{[ar]) A X T 
as a — > 00, where X T is a 7-stable Levy process known 
as the parent process. If 7 = 2, the parent process is 
the classical Brownian motion. Both, the process U T and 
the process X T are indexed by random operational (in- 
ternal) time r. In order to find a particle position at 
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the observable time t, we have to introduce the notion 
of the inverse-time a-stable subordinator St relating the 
internal and the observable times 



a- a N at 4 S t = inf{r : U T > t} 



(2) 



as a — > oo. Then, as a — > oo, the continuous limit of the 
CTRW process (UJ obtains the following form 

a- a ^R(N at ) a (a a )-^ R([a a S t ]) 4 X(S t ), (3) 

known as the anomalous diffusion process (llj . directed 
by the inverse a-stable subordinator St. It should be 
pointed out that the process U T does not have any fi- 
nite p-moments for p < a. Therefore, the subdiffusion 
is characterized by a power mean square displacement in 
time [lEEEl- 

However, there are physical phenomena, for example, 
the random motion of BPs in intergranular lanes on the 
Sun, where it would be desirable to get also a model that 
overcomes the infinite-moment difficulty while preserv- 
ing the subdiffusive behavior for short times [13j. The 
remedy was first proposed in the physical literature by 
Mantegna and Stanley 1J] . Their idea of truncated Levy 
flights served as a model for random phenomena which 
exhibit at small scales properties similar to those of Levy 
flights, but have distributions which at large scales have 
cutoffs and thus have finite moments of any order. Kopo- 
nen [15J , building on Mantegna and Stanley's ideas, de- 
fined the smoothly truncated Levy flights which stressed 
the advantage of a nice analytic form. Independently, 
the same family of distributions was described earlier 
by Hougaard [16] in the context of a biological appli- 
cation. However, different methods for the truncation 
were suggested also in the economic and statistical sci- 
ences jl7H19l |. but until the Rosihski's paper [H[ there 
was a lack of invariance under linear transformations for 
the distributions introduced, a significant property that 
the a-stable laws possess. He succeeded in finding the 
appropriate class of tempered stable distributions and 
processes @, Ell- 
in what follows we discuss properties of a diffusion 
process which is related to an inverse tempered a-stable 
subordinator. The Laplace image for the probability den- 
sity function (p.d.f.) of a tempered non-negative a-stable 
variable is 



f(u) = exp(6 a -(u + S) a ) , 



(4) 



where S is a positive constant and < a < 1 [6]. If S 
equals zero, the tempered a-stable p.d.f. becomes simply 
a-stable. Eq.(j4|) describes probabilistic properties of the 
tempered a-stable Levy process U(t) which generalizes 
the above mentioned process U T . 

Next, we will find its inverse process S(t) as in ©, 
where U(t) substitutes U T . If f(t,r) is the p.d.f. of 



U(t), then the p.d.f. g(r,t) of its inverse S(i) can be 
represented as 

9(r,t) = -^ f f(t',T)dt'. 

Taking the Laplace transform of g(r, t) with respect to t, 
we get 



(5) 



When u > 1 (t < 1) or S -> 0, Eq. © tends to 

g(r,u)=u a - 1 e-™ a , (6) 

which is the Laplace image of an inverse a-stable p.d.f. 
typical for subdiffusion. If u <C 1 (t 3> 1) or a — > 1, then 
Eq. ([5|) becomes the Laplace image of the Dirac delta- 
function. It follows from Eq. ([6]) that the p.d.f. of the 
inverse a-stable process is 



9(r,t) = — 

Z7TI 



J\lt—TU n ,OC— 1 



u a - L du = t- a F a (T/t a ), 



Br 



where Br denotes the Bromwich path, and the function 
F a (z) is a specific case of the Wright function (9L |20|. 

Subordination by an inverse tempered a-stable 
process. — Let the parent process X(t) have the p.d.f. 
h(x,r). Then the p.d.f. of the subordinated process 
X[S(t)] obeys the integral relationship between the pro- 
bability densities of the parent and directing processes, 
X(t) and S(t), respectively, 



p(x,t) = / h(x,r)g(T,t)dr. 



(7) 



In the Laplace space the probability density p(x, t) has 
the most simple form. Taking into account Eq. ([5]). the 
Laplace transform of Eq. (0 with respect to t gives 

p(x,u) = ±—^ h(x,(u + 5) a -S a ). (8) 

u 

For 8 = the latter expression becomes u a ~ 1 h(x,u a ). 

It is not difficult to calculate the moments of the pro- 
cess X[S(t)] if the moments of the process X(t) are 
known. For example, for the Gaussian process (7 = 2) 
the second moment is (X 2 {t)) = Dt, where D is a dif- 
fusive constant. Then the mean square displacement of 
X[S(t)] can be written as 

/•OO 

(X 2 [S(t)})= / (X 2 (r)) g(r,t)dr. 
Jo 

The Laplace image (XI) of (X 2 [S(t)}) has the form 

D 



{XI) 



u[(u + S) a - S a ] 



(9) 
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FIG. 1: (Color online) Mean square displacement of anoma- 
lous diffusion subordinated by an inverse tempered a-stable 
process. 



Consequently, the inverse Laplace transform of Eq. 
reads 



C- 1 



U[(« + £)<*- *-j 

where 



- Sy y a - l E a , a {8 a y a )dy, 



(10) 



k=Q 



T{ak + f3)' 



a > 0, (3 > 



is the Mittag-Leffler function [20|. The function (|T0|) 
gives rise to interesting asymptotic properties of the mean 
square displacement (X 2 [S(t)]). For t<l this displace- 
ment behaves as Di Q /r(cv + l), but for £ ^> 1 it increases 
linearly in time Dt/a (see Fig. [1]). Thus the anomalous 
diffusion, governed by the inverse tempered a-stable sub- 
ordinator, occupies an intermediate place between sub- 
diffusion and the normal diffusion. For short times it be- 
haves as subdiffusion whereas for the long times it resem- 
bles the properties of the normal diffusion. Let us call the 
diffusion subordinated by the inverse tempered a-stable 
process as a "tempered subdiffusion" . As is well known 
[2l| , the inverse a-stable process accounts for the amount 
of time, when a walker does not participate in a motion. 
By analogy, we may conclude that the process S(t) for the 
tempered subdiffusion represents a case, when a walker 
does not participate in a motion only for restricted in- 
tervals of time. At large time scales the walker begins to 
move randomly without any stopping as if a = 1. 

Equation of tempered subdiffusion. — Let L(x) be a 
time- independent Fokkcr-Planck operator (FPE), whose 
exact form is not important here. Let the ordinary FPE 
dh(x,T)/dr = L(x)h(x,r) describe the evolution of a 
particle subject to the operation time r. Acting by the 



operator L(x) on the image p(x, u) from Eq. (|8j). we find 



L(x)p(u,x) = [{u + S) a - S a ]p(x,u) 

- goo 



(11) 



where q(x) is an initial condition. When S = 0, the 
inverse Laplace transform of the latter expression gives a 
fractional representation of the FPE [l|, 0] 

p(x,t) = q(x) + 



dT{t- T ) a - 1 L{x) P {x,t). (12) 



r O) Jo 



In the case of the tempered subdiffusion the kernel in the 
integral representation of the FPE will be more complex, 
containing as a special case the kernel of Eq. (fT2"| for 5 —> 
0. Using the formal integral representation of the FPE 



p(x,t)=q(x)+ drM(t-T)L(x)p(x,T). (13) 
J o 

and taking the inverse Laplace transform of Eq. (fTTj) , we 
obtain the explicit form of the kernel M(t), namely 

M(t) = e- 5t t a - 1 E a , a (5 a t a ). (14) 

For t <C 1 (or 6 — > 0) this function takes the power form 
t a /T(a) as the kernel in Eq. (fl"2]) . However, for t 3> 1 (or 
a — > 1) M(t) becomes constant and, as a result, Eq. (fT3|) 
transforms into the integral form of the ordinary FPE. 

Tempered relaxation. — The commonly accepted theo- 
retical approaches to model relaxation phenomena as- 
sume [l| decay of an excitation undergoing diffusion in 
the system under consideration. In this framework, the 
relaxation function <p(t) describes the temporal decay of 
a given mode k and can be expressed [l0( through the 
Fourier transform of the diffusion process X[S(t)] 



4>{t) 



= ( e -kX[S{t)] 



(15) 



Here k > has the physical meaning of a wave number 
(the Fourier image of spatial coordinates). Starting with 
Eq. ([5]), we can write the Laplace image of Eq. (fT5")) as 



[{u + S) a -8 a ] 



(16) 



zt[$(fc) + (u + S) a - S a ] 

where $(fc) is the logarithm of the characteristic function 
of the process X(t). 

To expose the characteristic properties of the "tem- 
pered relaxation" we use the frequency-domain descrip- 
tion 



22. 23] of the relaxation phenonenon 



-*m dt. 

dt J 



(17) 



Then, for the relaxation under the inverse tempered te- 
stable process the function (fTT)) takes the form 



1 



a a + (ioj/uj p + <r) a ' 



(18) 



4 




KT 3 icr z io _1 io° io' io 2 10 3 



<0 



portant that they both can be considered on the unique 
base following the theory of subordinated random pro- 
cesses. We have derived a tempered form of the FPE 
and the relaxation function, as well as we have calculated 
the mean square displacement describing the processes. 
In Fig. [3J as an example, the propagator p(x, t) for the 
tempered diffusion with a — 2/3 and 5 — 0.5, is drawn. 
The cusp shape of the p.d.f. disappears when time in- 
creases. Thus our model occupies an intermediate place 
between subdiffusion and normal diffusion. We expect 
that our results will yield insights into the coexistence of 
subdiffusion and normal diffusion in nature. 

AS is grateful to the Institute of Physics, Wroclaw Uni- 
versity of Technology for kind hospitality during his visit. 



FIG. 2: (Color online) Absorption (imaginary) term of the 
frequency-domain relaxation function x( u ) = X { u ) ~ *X ( w ) 
for a — 0.75 . 
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FIG. 3: (Color online) Propagator p(x,t) for the tempered 
diffusion with a — 2/3 and S = 0.5, drawn for consecutive 
dimensionless instances of time t = 0.1, 1, 3, 10. The cusp 
shape of the p.d.f. disappears. 



where < a < 1 is a constant, and w p is the characteristic 
frequency of the relaxing system. 

According to Eq. (fl8|) . for a = the relaxation func- 
tion describes the Cole-Cole law. If a = I, the func- 
tion becomes exponential. In the case of a = 1 it has 
the Cole-Davidson form. The relaxation directed by the 
inverse tempered a-stable process takes an intermediate 
place between the superslow relaxation and the exponen- 
tial one (see Fig. [2j. Such a type of evolution is observed 
in relaxation experiments (see, for example, [22|). 

Conclusions. — In summary we have developed a novel 
approach to anomalous diffusion and nonexponential re- 
laxation from tempered a-stable processes. The model is 
broader than the purely subdiffusive case. It is very im- 
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